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I. INTRODUCTION 



Scattering by gravitational sources has been an important test for General Relativity since 
its inception. An early success for Einstein's theory came with verification of the prediction 
that starlight is deflected as it passes close to the Sun. The experiments conducted by 
Eddington and his team during the solar eclipse of 1919 were important in establishing GR 
as a credible theory of gravitation. More recently, gravitational lensing observations provide 
a strong test of the theory, as well as information about lensing mass distributions and 
distant astronomical sources. 

Gravitational scattering from astrophysical objects may be understood entirely in terms 
of geometric optics through the analysis of classical geodesies. Nevertheless, to obtain a 
deeper theoretical understanding of extreme objects such as black holes many authors have 
also considered the scattering of coherent waves. If the wavelength of the incident wave is 
comparable with the size of the event horizon then the wave will be diffracted by the black 
hole. Diffraction effects are responsible for many interesting phenomena in nature, such as 
glories and rainbows, so wave scattering from black holes is an interesting field in its own 
right, even if observations may not be realisable in practice. 

In this paper, we study how quantum waves are scattered by a Schwarzschild black hole. 
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17|. Most 



This problem has been considered in numerous papers and books 

authors study the massless scalar wave (with spin s = 0) for its mathematical simplicity, or 
massless electromagnetic (s = 1) and gravitational (s = 2) waves for physical relevance. In 
this paper, we focus instead on fermion scattering (s = 1/2), which has been less frequently 
discussed |l6j. We investigate the effect of particle mass on scattering cross sections and 
polarization. 

The interaction of an incident plane wave with a black hole may be understood in terms 
of three quantities: the total absorption cross section a^, the differential scattering cross 
section da/dQ, and the polarization V that is induced in initially unpolarized beam (note 
that ^ and V are functions of the scattering angle 6). In this paper we calculate these 
quantities using a variety of analytical and numerical methods. 

The paper is organised as follows. In section |n] we review classical scattering on the 
Schwarzschild background. We derive a new approximation for near-horizon deflection of 
timelike geodesies, and investigate classical spin precession on backward-scattering orbits. 
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In section ITTTl we describe analytic methods that can be apphed to fermion scattering. The 
perturbation-theory approach introduced by Doran and Lasenby is taken to the next 
order, which provides a new estimate of the cross section and polarization at low couplings. 
The glory-scattering approximation of Zhang and DeWitt-Morette Q] is extended to the 
massive case, employing the classical results of section IHl In section IIVI we discuss partial 
wave scattering theory, formulate the Dirac equation, and outline a simple numerical method 
for calculating phase shifts. In section|3we present the results of our numerical calculations, 
and compare with theory. We conclude in section IVII by discussing the significance of our 
results and possible future work. 

Note that in this article we employ units in which the speed of light (c), the gravitational 
constant (G), and Planck's constant over 2it (h) are set to unity. 

II. CLASSICAL SCATTERING 

We begin by reviewing some results for classical scattering and absorption on a 
Schwarzschild background. We then investigate the precession of classical spin vectors along 
scattering trajectories. As well as being of interest in their own right, the results of this 
section will aid our interpretation of the wave scattering cross sections presented in section 
IVl The results of this section are also required in section UTTl to derive a new approximation 
for glory scattering of massive spinor waves. 

The Schwarzschild spacetime is invariant under time-displacement and rotation. These 
symmetries imply the existence of Killing vectors, and the quantities conjugate to those 
Killing vectors are conserved along geodesies. That is, the particle energy E, and angular 
momentum L, are constants of motion. 

The motion of a particle of rest mass m on a Schwarzschild background is described by 
the orbital equation 



where u = 1/r, and is the polar angle. For descriptive purposes it is useful to define three 
further constants: the momentum p = [E'^ — m^)^/^, the impact parameter h = L/p, and the 
speed V = p/E. At large distances from the hole, each has a simple physical interpretation: 
p is the incident particle momentum, b is the orthogonal distance to the scattering centre. 




(1) 
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and V is the particle speed in units of c. Note that for null geodesies, v = 1 and p = E. 
With these definitions we may rewrite the orbital equation as 

fduV o 2 2M(1 1 

The unbound scattering solutions of (0) can be written in terms of elliptic integrals. For a 
comprehensive treatment, see Chandrasekhar |l5| . 



A. Deflection-angle Approximations 

In the weak-field limit, r ^ 2M, the deflection angle for timelike geodesies is approxi- 
mately 



For null geodesies we set f = 1 to recover the Einstein deflection angle for the bending of 
light, A0 ^ 4M/b. In the strong-field regime, the above approximation underestimates the 
deflection; incoming geodesies may be scattered through large angles or even orbit the hole 
several times before escaping (see Fig. Q). Trajectories that make too close an approach to 
the hole spiral inwards, and end on the singularity. The last geodesic to avoid the singularity 
defines a critical impact parameter be- All geodesies with b > be are scattered, whereas all 

^desics with b < be are absorbed. By considering the zeros of 0, one can show (see e.g. 

isl) that 

be = i^v^ + 2Qv^ -! + (! + 8i;2)3/2) ^'^ . (4) 

For null geodesies, the critical impact parameter is be = 3a/3M. The classical absorption 
cross section is simply the area of the circle defined by the critical impact parameter, a a = 

Tlbe". 

Many years ago, Darwin derived an approximation for the deflection of null (massless) 
rays that pass close to the critical orbit. He found that the deflection angle G is related to 
the impact parameter b by the equation 

& - &e = 216 X 3 V3M I ^ ~ M e-"e-® ^ 3.48Me-®. (5) 

Using the same techniques, an approximation for the deflection of timelike (massive) 
geodesies can be derived. This calculation is outlined in Appendix A, and the result is 




FIG. 1: (a) Photon geodesies around a Schwarzschild black hole, M = 1, close to the critical impact 
parameter be = SVSM 5.196M. The inner circle shows the event horizon at r = 2M, and the 
outer circle shows the unstable photon orbit at r = 3M. (b) Deflection angle i;^ as a function of 
impact parameter b. The plot compares the scattering angle with (i) the Einstein deflection angle 
4M/b, valid when b ^ be, and (ii) Darwin's approximation, b w SVSM + 3.48Me~'^, valid when 
b ~ be- 



b-be = M 
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(6) 



where Cc = (1 + Si;^)!/^. This result reduces to equation (0) in the case v = 1. 

A classical scattering cross section may be derived by considering a stream of parallel 
incident geodesies, as shown by Collins, Delbourgo and Williams 20|. Geodesies passing 
close to the critical orbit are scattered through angles that may be many multiples of 27r, so 
the total cross section is an infinite sum, 

db 



da 

dn 



sin 9 



d9 



(7) 



where the sum is taken over angles 6, 2tt — 6,2tt + 6, An — 6, etc. The cross section is 
divergent in the forward direction, da/dQ ^ 4f~^(l + v'^)^M'^6~'^ |2fl], due to the long-range 
nature of the interaction. It is also divergent (but integrable) in the backward direction, 
da/dQ oc Itt — 6\^^, because the element of solid angle tends to zero on-axis. The behaviour 
of the cross section close to 6* = and 9 = n may be understood using the approximations Q 
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FIG. 2: Classical scattering cross sections for a range of particle speeds. The cross sections diverge 
on-axis at 6^ = and 6* = vr as ^ Av'^^il + and ^ ~ k - e^]"^ respectively. 

and © given above, and at intermediate angles it may be found by numerically integrating 
the elliptic equations. The result is plotted in Fig. |21for a range of particle speeds. 

B. Spin Precession 

In the next section pilll we derive a semi-classical approximation to the massive fermion 
scattering cross section in the backward direction. To do this, we must first consider how 
the spin of the particle is affected by the gravitational interaction. For a massless wave, 
this is straightforward, as the helicity of the Dirac spinor remains unchanged. For a massive 
wave, the situation is substantially more complicated. However, progress can be made by 
considering the precession of the classical spin vector. We can then apply a semi-classical 
argument to relate the classical spin vector to the Dirac spinor itself in the A -C limit. 

We begin by introducing an angular momentum four-vector with components 
[s^ , , , s^], where t,r,6,(j) are the coordinates of the standard Schwarzschild metric Qf^^. 
The four- vector is spacelike, normalised, and orthogonal to the velocity v'^. Classically, 
angular momentum is parallel-transported along the geodesic, that is, s'^ + T^i.xv'^s'^ = 0. 
Inserting the standard Schwarzschild connection coefficients, and applying the orthogonality 
condition, we find the precession equations can be written as a pair of coupled first order 
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equations 

where s'^ = rs'^. Alternatively they can be combined in a simple second-order equation 



Now let us investigate a specific case, in which the spin vector is initially aligned with 
the direction of motion. At infinity, where u = 0, we start with initial conditions s"^ = —7 
and s'^ = 0, where 7 = (1 — v'^)^^^'^. The solutions to the precession equations may then be 
written as 

2 /„2\l/2, 



s 



(m) = — 7 (1 + ^^^cosx(^i) 



- {u/t]^) (1 + u^T]^) ^''^ {2Mu^ -u^ + 2Mr]\ + ^''^ sin (10) 

where i] = {'~fbv)~^ and 

X{u) = 7 / (1 + w^T]'^) {2Mw^ -w^ + 2Mrfw + 1/6^) '^^^ dw. (11) 
Jo 

For comparison, the scattering angle can be written in similar form as 



0(n) 



/ {2Mw^ -w^ + 2Mr]^w + l/b'^y^'^ dw. (12) 
Jo 



Equations (fTT| and (fT2|l can be expressed in terms of elliptic functions, but we prefer to 
solve them numerically here. 

In the next section we discuss glory scattering, which is related to the classical orbits 
that emerge in the backwards direction. We are particularly interested in the angle ^ that 
the classical spin vector makes with the outgoing radial direction, in the rest frame of the 
particle. For a clockwise orbit, we define 



where Xoo = xi^o), and mq = 1/ro is found from the radius of closest approach for the 180° 
orbit. The vectors and angles in Eq. (fT!?|l are illustrated in Fig. Efa). Figure Efb) plots ^ 
for the backward-scattering orbit as a function of particle velocity, and demonstrates that ^ 
decreases monotonically as v increases, from ,^ = 7ratf = 0, to,^ = Oatw = l. 
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FIG. 3: Spin precession around a 180° orbit, (a) Depicting the precession angle ^. Here, pi and p/ 
are the initial and final momentum directions, and Sj and sj are the initial and final spin directions 
in the rest frame of the particle, (b) The precession angle (solid line) as a function of particle 
velocity v. The dotted line is for reference only. 

III. SPINOR SCATTERING: ANALYTIC RESULTS 

In this section we review some approaches to calculating the fermionic cross sections for 
black hole scattering. The appropriate method depends on the ratio of the black hole size 
to particle wavelength. A convenient dimensionless measure of the gravitational coupling is 

-^^^ (-) 

where is the Schwarzschild radius, and X = h/p is the wavelength of the quantum particle. 
Reverting to natural units, we see that e = ME, which we use to label our figures. 

A. Low-coupling scattering cross sections: Perturbation theory 

In a recent paper Q|, Doran and Lasenby showed that the fermion scattering cross 
section can be expanded as a perturbation series. In effect, they calculated the gravitational 
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analogue of the Mott formula for scattering from a Coulomb potential. The series takes the 
form ^ 

^ = (^) [«o(^, e) + ea^iv, 9) + e'a,{v, 9) + ...], (15) 

with e as defined above, and ai{y,9) dimensionless functions. The perturbation series ap- 
proach is most appropriate in the long- wavelength limit, when e <^ 1. Doran and Lasenby 
showed that the first-order term oq is gauge-invariant and equal to 

1 + 2t;2 - 3t;2 s\n^{9/2) + - sm^{9/2) 
Av^ sin^(^/2) 

In the massless limit {v 1), this reduces to the unpolarized Mott scattering formula. The 
cross section will tend to the first-order result in the long wavelength limit, as r5/A 0. 
In Appendix B we show how to calculate the next term in the perturbation series in the 
Kerr-Schild gauge. We find 

TT [(3 + Av^ + v%2 - I sin(g/2)|) - {Iv^ + v') sin^(g/2)] 
^) = 4.3|sin(g/2)|3 • (1^) 



, , ± -r — ou bill <yuiz,)^ v —V &m yu/^) . 

ao[V,U) = .../t„: 4//i/o\ ■ y^^J 



The equivalent result for Coulomb scattering was successfully calculated by Dalitz |2J| many 
years ago. The second order matrix element also allows the calculation of the degree of 
polarization V induced in an initially unpolarized beam. The net polarization is in the 
direction orthogonal to the plane of scattering, Pj x pj. To first order, 

^ fGMm\ v{3 + v')t^n'{9/2) ln(sin^(g/2)) . 

V he J 2[l + 2v^-3v^sm^{9/2)+v^-vHm\9/2)]'^''' ^ ' 

where we have included the dimensional constants for clarity. The confirmation of the 
gauge-invariance of results ()17p and p8|) awaits further work. In section |3 we compare the 
perturbation series and polarization with numerical results. 

B. Higher-coupling scattering cross sections: Glory and spiral scattering 



Many authors 8|, |lO|, |ll| have noted that black hole wave scattering produces diffraction 
effects that are familiar from optical phenomena. Two such effects are glory and spiral 
scattering. A glory is a bright spot or halo that appears on-axis in the backward direction 
from the scatterer. Spiral scattering, or orbiting, creates oscillations in scattered intensity 
at intermediate scattering angles. In a general analysis of semi-classical scattering. Ford and 



Wheeler 



22| showed that both effects may be understood with reference to the deflection of 



classical paths. A glory occurs whenever the deflection angle passes through a multiple of 
TT, and spiral scattering occurs when paths orbit the scattering centre multiple times. For a 
discussion of the application of semi-classical techniques to the black hole case, see lll| or 

Q. 



Ford and Wheeler [2^ derived a semi-classical approximation of the glory scattering cross 
section for scalar (s = 0) waves. This approximation was extended to arbitrary spins by 
Zhang and DeWitt-Morette in |9|] using path integral methods. For massless waves, the 
backward glory cross section is approximated by 

db 



da 

dn 



d9 



J2s\EbgSme), (19) 

glory 

where s is the spin of the particle, is a Bessel function, and bg is the impact parameter 
at which 6 = tt. In combination with Darwin's result (0) for the deflection of a massless 
particle that passes close to the horizon, we have bg ^ 5.3465M, and \db/d6\ = 0.1504M, so 



, ^ da 

dn 



^ 2ttEM X 4.30 X J2^,(5.3465EMsin^). (20) 

glory 

For the scalar wave, the intensity has a peak in the backwards direction, whereas for the 



spinor wave the intensity is zero on-axis. 

We now derive a similar approximation for the massive spin-half case, by taking into 
account the effects of classical spin precession. Let us begin by considering a backward- 
scattering orbit deflned by pj = i = — p/, that lies in the (pz plane. As we saw in Section 
111 Bl the classical spin vector is rotated by an angle of 2xoo around this orbit. The quantum 
spin vector, on the other hand, is found by the double-sided action of the Dirac spinor, 

= ('?M7'^75|'?A), where 7^ are the standard Dirac-Pauli matrices. Assuming that the 
classical and quantum spin vectors coincide in the classical limit, A/r^ 0, we conclude 
that, in the rest frame, the Dirac spinor must transform according to 

tpf oc (^cos(xoo)/ + sin(xoo)7<^73) 4'i (21) 

where J is the identity matrix, and 7,/, = — sin07i -|- cos 072. Now consider the cross section 
on-axis. Following the arguments of Zhang and DeWitt-Morette ^9|, we note that there is a 
circular degeneracy here, and that paths from all scattering planes will interfere coherently. 
Around the circle < < 27r, the 7^73 term of (j?H) sums to zero, and the intensity on-axis is 
therefore proportional to cos^(xoo)- Since it is the spinors, rather than the vectors, that are 
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summed coherently, the half-angle of precession arises naturally. The arguments of Zhang 
and DeWitt-Morette can be extended in a straightforward way to derive an approximation 
for angles close to tt, leading to 

cos^iXoo)Jo\Ebvsm9) ^^^^ 
+ sm'^{xoc)Ji^iEbvsme) 

Finally, we may use (jH)) to evaluate 6 at = tt, and substitute into the above result. 
Alternatively, we may choose to evaluate bg numerically for a more accurate result. In 
section IVl we compare both approaches with our numerical results. 



glory d9 



IV. SPINOR WAVE SCATTERING 



Before considering the Dirac equation on a Schwarzschild black hole background, we first 
discuss the general theory of relativistic spin-half scattering from a spherically-symmetric 
source. 



A. Scattering from a spherically-symmetric source 



In considering spherically-symmetric scattering we will employ the two-component spher- 
ical spinors Xk that are eigenvectors of the angular equation 



{cT-L + l)x^, = Kxl 



(23) 



where L = [L^, Ly, Lz] is the angular momentum operator, and the components of cr = 
[ax, (Ty, az] are Pauli spin matrices. Angular states are labelled by the eigenvalue k which is 
related to the overall angular momentum j by 



« = ±(J + 2^ 



/ + 1 where / = j - 1/2 
— / where / = j + 1/2 



Note that our sign convention for k may differ from the literature 
spherical spinors are 



m=±l/2 



(24) 

Explicitly, the 
(25) 
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where 5'^ = C{jij2jtot', mim2) are the Clebsch-Gordan coefficients, ¥[^{6, 0) are spher- 

ical harmonics, and ^'^^ the spin-up and spin-down spinors, = [q] aiid = [i]- 
The positive and negative n spinors are related by 

^rX'. = X'-., (26) 

where = r^^ xVj. 

Consider, first, the free-particle Dirac equation, ij'^d^'^ — m\E' = 0. It admits plane wave 
solutions; for example, a plane wave propagating in the z direction can be written as 

*Sane =1 ,1 e^^'e-^''^ (27) 



E+m"^X 



It also admits separable solutions in spherical coordinates of the form 

^„_l/"i(r)x!:\ (28) 

^ \U2{r)X-^j 

where ui{r) and U2{r) are radial functions, and Xk the Pauli two-component spherical spinors 
()25|). The free particle radial equations can be written in matrix form as 

/t/r i{E + m) 
i{E — m) —n/r 

and their solutions are spherical Bessel functions 



d 1 






dr \ 




H 



(29) 



ur{r) = rji{pr), u^2 (r) = " rjj(pr), (30) 

£/ + 777, 

where I = Z + S*,;. The plane wave solution ()27|1 can be written as a sum of spherical solutions 
P^. as 

^^ane = 6^^* (4^) ^ (2/ + 1) |j ; Om) ^ ; '^'^ . (31) 
The asymptotic behaviour of the radial solutions is found by noting 

ji{pr) — sin ( pr J , as r ^ oo. (32) 



pr \ 2 

When an interaction is introduced, the unperturbed plane wave is no longer an eigenstate 
of the Hamiltonian. Instead we look for a time-independent solution with the asymptotic 
behaviour 

^ ^ ^piane + ^e^^^e"^^*, (33) 
12 



where C is a four-component spinor to be determined. Asymptotically, the wave function 
is interpreted as the sum of an incident plane wave and a radially-outgoing scattered wave. 
The effect of an interaction is to introduce phase shifts into the asymptotic solution (|32j). We 
define an interaction phase shift 5^ by comparing the asymptotic solution with the free-space 
form (jSg), 

u{'\r) — > sin ( pr — tt + ^n] , as r — >• cx3. (34) 



2 

An expansion of the form ()33|) is only strictly valid if the interaction is localised. In a long- 
range potential, such as the Coulomb 1/r potential, the asymptotic solutions are modified by 
the presence of a radially-dependent phase factor. The same is true in a gravitational field, 
as we show in the next section. The extra phase factor means that incident plane waves are 
distorted by the presence of the potential even at infinity. However, the extra phase shift 
is independent of k, so does not contribute to 5^ and the phase difference between partial 
waves. 

Let us now assume that the wave is incident along the z-axis, and is in a superposition 
of spin-up and spin-down states, with upper components 

E = f ""M • (35) 

We now construct a partial wave series with asymptotic form ()33p . For clarity we need only 
consider the upper components of the scattered wave, given by 

Upper(C) = J]c„5rx^ (36) 

mr 

where 

B7 = E(2^ + lY^'ie'^'^ - l)C(/lj; Om)C(/ij; m - r, r)r---(p^). (37) 

K 

1/2 

Using the properties of Clebsch-Gordan coefficients, it is straightforward to show that B-^^^ = 
B_y2- The scattered wave spinor can be written as Aj^m'^mX"^^ where A is a transition 
amplitude matrix, with components 

■4= I 5? ""['ig-zm+isw^-fl. (38) 

-° -1/2 -°-l/2 



13 



Here, / and g are complex amplitudes, and it can be shown (see, for example, 23[) that 
the rotation vector n is orthogonal to both the incident and scattered direction, fi = x 
Pf/\Pi ^ Pf\- The scattered intensity from an unpolarized beam is the sum of the squares 
of the amplitudes, 

(39) 

In terms of phase shifts, the scattering amplitudes are 



1=0 



^ 1=0 

where x = cos 6. There are two limiting cases to consider. First, in the non-relativistic limit 
we expect the interaction to depend only on orbital angular momentum (/). In this case, 
we expect 5^+1 ~ and the spin is then unchanged by the interaction, and we recover the 
non-relativistic scattering series {2ip)~^ J2i=oi'^^ + ^)e'^^^^Pi{x) (where r]i = Si+i = S-i). 

Second, in the highly-relativistic regime E ^ m, the phase shifts instead depend only on 
total angular momentum (j), and so 6^ ~ 6-^. Then, for 6^ 7^ 0, 

/(^) ^ i + 1)^'^'"^ [^'+i(^) + ^'(^)] (42) 

^ 1=0 

9(0) ^ ^ E ^'^'"^ [^^'(^) - ^ii(^)] (43) 
1=0 

Jsing the properties of the Legendre polynomials, it is straightforward to show (eg. Mott 



2J) that 

g = tan(^/2)/ (44) 

in the relativistic limit. Physically, this implies that the helicity of a massless wave remains 
unchanged by the interaction. 

If the particle has mass then the scattered beam may become partially polarized by the 
interaction. The degree of polarization V induced in an initially unpolarized beam is 

^ = (45) 

\f\ + 191 

in the direction n. This effect is known as Mott polarization, and is well known for the 
Coulomb interaction of an electron with a nucleus. Conversely, if the incident beam is already 
partially polarized (with polarization vector P) then the wave is scattered asymmetrically, 
with asymmetry factor VP ■ n. 
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B. The Dirac equation on a Schwarzschild background 



We now consider the Dirac equation on a Schwarzschild spacetime. Before starting, we 
must choose an appropriate coordinate system. A disadvantage of standard Schwarzschild 
coordinates is that they are valid only in the exterior region, r > 2M. Geodesies do not cross 
the horizon in a finite Schwarzschild coordinate time, and all solutions to the wave equation 
are singular at the horizon. Here, we prefer to work with coordinate systems that remove 
the coordinate singularity at r = 2M. Two possibilities are advanced Eddington-Finkelstein 
coordinates, with metric 

rfs^ = (1 - 2M/r)df - {AM/r)dt dr - {1 + 2M/r)dr^ - r^dn^, (46) 

a.dPai,Uev.-Gul.t.a.dcoo..di.atesny,with.et* 



ds'^ = {l- 2M/r) dt^ - ^m/rdt dr - dr^ - r^dn\ (47) 

Both coordinate systems deal smoothly with the horizon, and cover regions I and III of the 
Penrose diagram of the fully extended Kruskal manifold j^] • For a discussion of gauge- 



invariance and the choice of coordinates, see 

The time coordinate of the Painleve-GuUstrand (PG) system has a simple physical in- 
terpretation. It corresponds to the proper time as measured by a observer in freefall who 
starts from rest at infinity. For this reason, we now formulate the Dirac equation on a black 
hole background described by PG coordinates. 

We let {70,71,72,73} denote the gamma matrices in the Dirac-Pauli representation, and 
introduce spherical polar coordinates (r, 6, 0). From these we define the unit polar matrices 

7r = sin6'(cos07i + sin072) + cos6'73 , 
761 = cos6'(cos07i + sin072) — sin6'73 , 

70 = - sin0 7i + COS0 72. (48) 
In terms of these we define four position-dependent matrices {gt, dr, de, g<i>} by 



I2M 

9t = lQ + \ l — 7r , 9e = rje , 

gr = 7r, g<p = r sm9-f^. (49) 
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These satisfy the anti-commutation relations 



{9f^,9u}+ = 2g^ul (50) 

where Qf^^ is the Painleve-Gullstrand metric of equation (07j). The reciprocal matrices 
g^, g^, g^, g'^ are defined by the equation 

{g'^,g,U = 2Si^I, (51) 

and both sets are well-defined everywhere except at the origin. 
The Dirac equation for a spin-half particle of mass m is 

ig'^'^f.i^ = m^P, (52) 

where 

V;.^ = {d, + '-rfj:^p)tP, = ^[7„,7/3]- (53) 

The components of the spin connection are found in the standard way j3] and are particu- 
larly simple in the Painleve-Gullstrand gauge 



An advantage of our choice of metric is that the Dirac equation can now be written in a 
manifestly Hamiltonian form 

—) + i7 j ^ = (55) 

where ^ is the Dirac operator in flat Minkowski spacetime. The interaction term is non- 
Hermitian, as the singularity acts as a sink for probability density, making absorption pos- 
sible. 

The Dirac equation is separable in time, with solutions that go as exp{—iEt). The energy 
E conjugate to time-translation is independent of the chosen coordinate system, and has a 
physical definition in terms of the Killing time j27|. We exploit the spherical symmetry to 
separate the spinor into 

^^e-^fu.ir)M (56) 
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The trial function ()56|) results in a pair of coupled first-order equations 
d (uA I 1 v^2M/r 



(l-2M/r) , 
dr 





H 






(E- 


- m) - 



^2M/r 1 

K/r i{E + m) - {2M/ 

(2M/r)i/2(4r)"i -/t/r 



(57) 



The equations have regular singular points at the origin and horizon, and an irregular 
singular point at r = oo. Analytic solutions to radial equations of this nature have been 
investigated by Leaver ,29;] and others [s^. Here, we prefer to use series solutions around 
the singular points as initial data for a numerical integration scheme. 



C. Series Solutions and Boundary Conditions 

As is clear from (jFfjl . there is a regular singular point in the coupled equations at the 
horizon, r = 2M. We look for series solutions 

U^or =h]={r- 2My Y^rA[r- 2M)" (58) 

\U2 j „.=0 

where s is the lowest power in the series, and a„, 6„ are coefficients to be determined. On 
substituting into (jFfj) and setting r = 2M we obtain an eigenvalue equation for s, which has 
solutions 

s = and s = -\ + 4iME. (59) 

The regular root s = ensures that we can construct solutions which are finite and con- 
tinuous at the horizon. Regular solutions automatically have an ingoing current at the 
horizon. The singular branch gives rise to discontinuous, unnormalisable solutions with 



an outgoing current at the horizon 



We therefore restrict attention to the regular, 



physically-admissable solutions. The eigenvector for the regular solution is 

(60) 



1 ao 1 




[boj 





K-2iM{E + m) + 1/4 
K + 2iM{E -m) - 1/4 

In order to expand about infinity we need to take care of the irregular singularity present 
there. There are two sets of solutions, t/i^*^*^ and which asymptotically resemble 
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outgoing and ingoing radial waves with additional radially-dependent phase factors. To 
lowest order, 



.p/ {E + m) 



-p/{E + m) 



where the phase factors (t)2{r) are given by 



M 



01 (r) = EVSMr, 02 (r) = — (m^ + 2p^) In(pr). (62) 

p 

The general solution as r ^ cxd is a superposition of the ingoing and outgoing waves, 

. U^or r^2M 

for each angular mode. The magnitudes of A™ and ^4°"* determine the amount of scattered 
and absorbed radiation present. The phase shift 5^ is given by their ratio, 

/lout 

^2^S. = (-1)'+^-^. (64) 



D. Numerical Method 



To calculate the scattering amplitude we must first determine the phase shifts. Various 



analytical methods for calculating phase shifts have been success: 



uU y a pplied to the black 



hole scattering problem, for waves of spin s = and s = 2. [lOl lUl, ll^ |32] . Here we calculate 
the phase shifts using a numerical method. It is a relatively straightforward task to inte- 
grate the radial equation (jKTjl and apply the ingoing boundary conditions to compute 
the phase shifts defined by (jMjl . We determine the ingoing and outgoing coefficients A'-^ and 
^4°"* by matching the horizon solution onto the infinity solution. In practice, this involves 
choosing two boundary points r^^^ and rmax, close to the horizon and infinity respectively. 
We start at r^i^i with series solution f/hor('"min)5 and integrate out to rmax where we com- 
pare the numerical solution with the asymptotic form AJf ?7c^'^'' (rmax) + A'^^U^^^\rj^s,^) . The 
phase shift calculated by this method is independent of our choice of rmin and rmax, provided 
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that they are sufficiently close to the boundary points. The phase shift should also be inde- 
pendent of the choice of coordinates. We repeated our analysis using Eddington-Finkelstein 
coordinates to check the numerical accuracy of our results. 

For the scalar wave in the large-/ limit 16| the phase shifts approach "Newtonian" values 
which arise in a strict 1/r potential, 

r(/ + 1-^/3) , . 

r(/ + 1 + ^/3)' ^^""^ 

where (3 = M{2E'^ — m?)/p = ME{1 + v'^)/v. Phase shifts of this form also occur in the 
non-relativistic Coulomb scattering problem (see, for instance, 2^), with j3c = Zam/p. It 
is well known that the problem of non-relativistic Coulomb scattering can be solved exactly 
in parabolic coordinates (see, for example, 3^). The pure-Newtonian phase shift series 
defined by 



^ 1=0 



T{l + l-z(3) ^ 

r(/ + 1 + 2/5) 



Plicose) (66) 



can be summed using the Coulomb result to 



It is more difficult to sum the series directly, because it is poorly convergent. This is 
related to the fact that an infinite number of Legendre polynomials are required to describe 
the divergence at 6' = 0. Many years ago, Yennie, Ravenhall and Wilson js^l outlined one 
possible way around this problem. Given a Legendre polynomial series 

fi9) = J2^fPiioose) (68) 

1=0 

that is divergent at ^ = 0, one may define the mth reduced series, 

(1 - cos9rf{9) = ^aJ^^P; (cos ^). (69) 

The reduced series is obviously less divergent at = 0, so one may hope that the re- 
duced series converges more quickly. Using the properties of the Legendre polynomials, it is 
straightforward to show that the new coefficients aj"*"^ are related to the old coefficients a] 
by the iterative formula 

(j+l) (i) ^ + 1 (i) ^ (i) ,r,n\ 

a) — a, ; a)' ; a) V (70) 

We have found that this to be an excellent method for summing the series numerically, and 
that two or three iterations are sufficient. 
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V. NUMERICAL RESULTS 



In this section we present the results of our numerical calculations, and compare with 
theoretical expectations outlined in sections [H] and 11111 

A. Absorption Cross Sections 

The existence of an ingoing current at the horizon implies that flux is absorbed by the 
black hole. The absorption cross section cta is defined as ratio of the absorbed and incident 
fluxes. For the spin-half wave, 



At low energies, absorption is dominated by the j = s mode. For massless waves, the 
low energy limits are: (i) a a = IQnM^ for s = 0, (ii) o"a = 27rM^ for s = 1/2, (iii) 
CTA = fniMEfM^ for s = 1, and (iv) (Ta = ^TriME^M^ for s = 2 (see |35i]). Note 
that the low-energy scalar and spinor cross sections are proportional to the black hole area, 
whereas the electromagnetic and gravitational cross sections tend to zero in this limit. 

The massless scalar and fermion cross sections are shown in the left panel of Fig |3] For 
massive particles, the cross section diverges in the low energy limit, as shown in the right 
panel of Fig 0] In all cases, the cross section approaches the geometrical optics value of 
o"A = 277rM^ in the high-energy limit. For more details on the absorption cross sections of 
massive spin-half particles, see Doran et al. j^^. 

B. Scattering at low couplings, A » 

In section IIIII we outlined our expectations for scattering at low couplings. We now 
test them numerically. Figure El compares the perturbation series fjlSj) with a numerically- 
calculated cross section, at ME = 0.04 and v = 0.6. It is clear that the approximation is 
improved by the inclusion of the second-order term, Eq. ()17|1 . 

Figure El compares the perturbation results with the numerical calculations across a range 
of values of ME and v. The first-order series fll(i|l is shown as a dotted line, and the 
second-order series ffT7|) is shown as a dashed line. As expected, we find that the second- 
order correction improves the fit, and the approximation is in excellent agreement when 




(71) 
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FIG. 4: Absorption cross section, (a) Massless. Shows the absorption cross section of a massless 
scalar [sohd] and a massless fermion [dashed]. The classical limit at 277r is also marked, (b) 
Massive. The solid line shows the quantum absorption cross section for GmM/hc = 0.2, and the 
dotted line is the classical cross section (jlj). 




180 



Scattering angle 9 (deg) 

FIG. 5: Comparison of the numerical cross section [solid] with the first-order [dotted] and second- 
order [dashed] Born approximation at ME = 0.04, v = 0.6. Note that the scale on the y-axis of 
the left plot is logarithmic. 

MEv ^0.1. At higher couplings, the truncated perturbation series is still accurate for 
small scattering angles, but it fails to predict glory and spiral scattering oscillations. In this 
regime, ~ A, a full numerical calculation is necessary. 

The perturbation calculation also predicts the magnitude of Mott polarization, which we 



21 




30 60 90 120 150 180 30 60 90 120 150 180 30 60 90 120 150 180 





30 60 90 120 150 180 



FIG. 6: Comparison of the numerical cross section [solid] with the first-order [dotted] and second- 
order [dashed] Born approximation across a range of ME and v. The x-axis shows the scattering 
angle (deg) and the y-axis shows the cross section, log^o (-^~^^)- 

now test against numerical results. The first-order polarization is given by ()18p. whereas the 
numerical polarization is determined by ()45|1 . Figure [3 compares the perturbation series and 
numerical results for a low coupling, ME = 0.002. The left-hand plot (a) shows that there is 
an order-of-magnitude agreement between the two results. Naively, in the limit ME 0, we 
would expect to find perfect agreement between the two. We do not observe this. However, 
it is remarkable that we can obtain excellent agreement by artificially removing the effects 
of absorption from the numerical results. If we set the imaginary part of the numerical 
phase shifts to zero, we obtain the right-hand plot. Fig. Efb). It is apparent therefore that 
the perturbation method is somewhat flawed, as it does not account for absorptive effects. 
This is related to the fact that, whereas plane waves form a complete basis for the hydrogen 
atom, they do not form a complete basis on a black-hole spacetime. We discuss this problem 
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FIG. 7: Comparison of the numerical polarization [solid] and the perturbation calculation (|18|1 
[dotted] for ME = 0.002. (a) With absorption. The perturbation result ((TH|) does not match the 
numerical result, (b) Without absorption. Enforcing the artifical condition that numerical phase 
shifts are purely real removes the effect of absorption (see text). We then observe an excellent 
agreement with the perturbation result. 

further in Section IVTl 



C. Scattering at higher couplings, A ~ 

Figure IHl shows numerically-determined scattering cross sections for a massless fermion at 
larger couplings. As expected, we observe spiral scattering oscillations, and a zero on-axis 
in the backward direction. The magnitude of the oscillations, and their angular frequency, 
increases with black hole mass. The zero in the backward direction can be verified by sub- 
wave series The massless glory approximation 
proves a good fit to the numerical results close to 
6 = 180°, as shown in Fig. El It correctly predicts the approximate magnitude and angular 
width of the oscillations. 

The scattering cross sections for massless scalar and spinor waves are compared in Fig. 
ITUl In the backward direction the scalar wave is at a maximum on-axis, whereas the spinor 
wave is at a minimum, as predicted by the glory approximation. At intermediate angles, 
the regular oscillations in the intensity of the scalar and spinor waves appear to be 180° out 
of phase. 



stituting Pii^n) = (—1)' into the partia 
(1201) of Zhang and DeWitt-Morette ^ 
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ME = 1.0 




-2 I 1 ^ 1 ^ 1 ^ ^ 1 1 

45 60 75 90 105 120 135 150 165 180 

angle 9 (deg) 

FIG. 8: Glory scattering cross sections for the massless spin-half wave for various couplings ME = 
GME/hc^. 

Scattering cross sections for the massive spin-half wave are shown in Fig. El Here, 
the solid line shows the partial- wave result, and the broken lines show the semi-classical 
approximations to the glory (j22j) . The dotted line is the result of using the modified Darwin 
approximation ©, and the dashed line is calculated by integrating the orbit equations 
numerically. It is clear that the numerical approach (dashed) provides the better fit. This 
is because, for slower particles, we find bgiory is too far from for the first-order Darwin 
approximation to be sufficient. In general, the Darwin approximation underestimates the 
magnitude of the glory. 

Figures IT^ a) and ll2r b) show the magnitude of Mott polarization as a function of scatter- 
ing angle, for a range of couplings. The relationship between polarization, scattering angle, 
particle velocity, and black hole mass is not immediately obvious. However, Fig. shows 
that the oscillations in polarization are clearly related to the glory and spiral scattering 
oscillations, suggesting a semiclassical model for Mott polarization may be possible. This 
awaits further work. 
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120 




scattering angle 6 

FIG. 9: Glory oscillations for the massless spinor and scalar waves. This plot compares the 
numerical glory oscillations (solid) with the approximate formula in the text (|2())1 (dotted), (a) 
Scalar, s = 0. (b) Spinor, s = 1/2. 
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angle (deg) 
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45 60 75 90 105 120 135 150 165 180 

angle 9 (deg) 

FIG. 10: Scattering cross sections for the massless scalar [solid] and spinor [dashed] waves at a 
coupling of ME = 4. Note that the bottom plot shows the same data as the top plot but with a 
logarithmic scale on the y-axis. 
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Scattering angle 9 (deg) 



FIG. 11: Glory oscillations for the massive spin-half wave, over a range of particle speeds v. 
The solid line shows the numerical cross section. The dashed line shows the result of the semi- 
classical approximation H22|) . and the dotted line shows the result of using the modified Darwin 
approximation to estimate the glory impact parameter. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we have presented an investigation into the scattering and absorption of 
massive spin-half waves by a Schwarzschild black hole. We extended the approximations 
of Darwin and Zhang and DeWitt-Morette ^ to derive a semiclassical (A <^ r^) ap- 
proximation for massive fermion glory scattering. We also took the perturbation analysis of 
Doran and Lasenby 14] to next order, to derive new formulae for the scattering cross section 
and Mott polarization in the low coupling limit, A ^ r^. Finally, we applied partial- wave 
techniques to solve the Dirac equation on the Schwarschild spacetime, and presented our 
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FIG. 12: Mott polarization as a function of scattering angle for a range of ME. Waves with mass 
are partially polarized by the interaction. As the coupling ME is increased, oscillations arise in 
the polarization, (a) For v = 0.9. (b) For v = 0.2. 
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FIG. 13: Glory oscillations and polarization for ME = 2 and v = 0.9. (a) Partial polarization, 
P. (b) Scattering cross section. The solid line shows the total cross section, and the dotted and 
dashed lines show the contributions from the amplitudes / and g defined by Eq. (|40() and (|41() . 

numerical results. We showed that our approximations provide good fits to the numerical 
results in the appropriate limits. 

We showed that non-zero particle mass leads to a phenomenon familiar from electromag- 
netic scattering: Mott polarization. This is a partial polarization created in the direction 
orthogonal to the scattering plane. In the electromagnetic case, an estimate of the polariza- 
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tion arises from the second-order term in the Born approximation [2l|. In the gravitational 
case, a similar estimate arises from our second-order perturbation analysis in the Kerr-Schild 
gauge, Eq. (jl8|) . However, our estimate only agrees with the numerical calculation if we 
artificially suppress black hole absorption (see Fig. [7j). This highlights a shortcoming in the 
Born approximation method: it cannot account for absorption effects. To see why, consider 
the Hamiltonian form of the Dirac equation in PG coordinates, Eq. ()55p. As r ^ 0, all 
solutions go as ~ ^-3/4^ ^^^^ have a net ingoing current. Clearly, such solutions cannot 
be described by a linear sum of plane waves used in the Born approximation. A related 
observation is that the Hamiltonian is not Hermitian as probability flux is removed at the 
origin It is not immediately clear how we extend the Born approximation to such open 
systems. Despite such difficulties, however, we have seen that the Born approximation in 
its present form provides a good fit to numerical scattering cross sections when A ^ r^, and 
it remains a useful tool. 

One intriguing question raised by this work is: when, if ever, are fermion diffraction 
patterns from black holes of physical relevance? Since diffraction effects are only significant 
when Ts ~ A, the question may be rephrased: do small black holes exist in nature? By 
small, we mean black holes with an event horizon smaller than the Compton wavelength of 
the lightest fermion (thought to be the electron neutrino). Recent experimental evidence 
suggests that neutrinos have non-zero mass. A neutrino with a hypothetical mass m^^ ~ 0.01 
eV has a Compton wavelength of 1 x 10"^ m, and this size corresponds to a black hole of 
mass 

M ~ 10^2 i^g^ Black holes of this magnitude or smaller are unlikely to be created by 
astrophysical processes; however, it is possible that small primordial black holes were created 
in the early universe. An alternative scenario for small black hole genesis arises in recently- 
proposed theories with Large Extra Dimensions If such theories are well-founded, it 
is possible that black holes may one day be created at a particle accelerator (see j3| for a 
review of these ideas). 

Despite these possibilities, fermion diffraction patterns may never be observed. A more 
physically-relevant problem is that of gravitational- wave scattering from rotating, large black 
holes (M > 10'^'^ kg). Previous studies of massless spin-two waves on the Kerr spacetime 
have shown that the situation is complicated by effects including frame-dragging, spin- 
rotation coupling, and superradiance. Futterman, Handler and Matzner present s = 2 
scattering cross sections for a special case: the alignment of black hole spin axis and incident 
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wave momentum. However, we are not aware of a general analysis of off-axis gravitational 



scattering, though progress has been made for the simpler spin-zero case j39|. We suggest 
that a general study of spin-half scattering from a Kerr black hole would be a useful addition 
to the literature, for two reasons. Firstly, Unruh has shown that superradiance does not 
occur for massless spin-half waves. Secondly, though we would expect to see the effects of 
spin-rotation coupling, the fewer degrees of freedom in the s = 1/2 case will make the results 
easier to interpret. Such a study would contribute to our understanding of gravitational- 
wave scattering, as well as being interesting in its own right. We therefore hope to extend 
the methods of this paper to the Kerr black hole in the near future. 



APPENDIX A: NEAR-HORIZON SCATTERING APPROXIMATION 

In Darwin derived an approximation for the deflection B of an unbound massless 
particle passing close to the horizon. Here we extend that calculation to the massive case, 
where f < 1. 

Following the approach of Chandrasekhar [l^, we look for roots of the orbit equation 
((21). For unbound orbits that avoid the singularity, the orbit equation has three three real 
roots, Ui < 0, M3 > U2 > 0. The roots can be written in terms of the eccentricity e and latus 
rectum / 

e-1 e + 1 1 2 ,^ , 

Ul = ; — , Uo = — ; — , Ms = — rr — T (All 

where e and / satisfy 

1 -/i(3 + e2) _l-v^ 
Ml f262 ' 

l(e^-l)(l-4/.) = l, (A2) 
and fi = M/l. Making the substitution u = {1 + ecosx)/l in (0) yields 

^ = (1 - + 2/ie) V2 [1 _ e cos'ixm] (A3) 

where 

e = ^-^ (A4) 

The variable x runs from < x ^ Xoo = cos^^(— e^^). The solutions of ()A3|) are elliptic 
functions, and we can write the final angle as 

2 



;i -6/i + 2/ie)i/2 
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[K{k)-F{l7r-lxo.,k)] (A5) 



The critical orbit occurs where U2 = M3, k = 1. Then 



= 2M(Vl + 8f2 + 3), 

^8^7^ + 20t;2 _ 1 + (1 + 8t;2)3/2M 



(A6) 



V2i 

The critical parameters are functions of the particle speed v only and lie in the ranges 
1 < Cc < 3, 8M < < 12M, and 3M < r^^^ < AM. We now consider orbits close to the 
critical orbit, by expanding in a power series in 6 = I — Ic 

I = lc + S, 

e = e. + 



b = bc + M 



e, + l)(e, + 3)M 



2y2t;(e, + 3)i/2(e, + 1)5/2 m2 



fc2 = l-— i -A. (AT) 

ee(e, + 3) M ^ ^ 

Note that 6 varies only quadratically with 6. We employ the following approximations for 
the elliptic functions 

K{k)^ In- ^ 



(1 - P)l/2 ' 

/ 1 + sin z 

2^^ 1 ■ — 

\ 1 — sm z 

and substitute into ()A5|) . Finally, we define the deflection 9 = 20oo — vr to obtain result 



F{zA)^Un(l±^) , (A8) 



APPENDIX B: SECOND-ORDER BORN SERIES 

Doran and Lasenby showed how to calculate the first-order contribution to the 
fermion cross section in a consistent way. By expanding the scattering amplitude as a 
power series in the black hole mass M, they found the first term in the scattering amplitude 
series to be 

Mi = Us{Pf)tiUr{Pi) , where f 1 = - ^^^^ (2^7° - m) (Bl) 

\Pf-Pi\ 
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and demonstrated its gauge invariance. Below, we calculate the second term in the scattering 
amplitude, A^2, in the Kerr-Schild gauge. We find the real part to be 



\Pf-Pi\ ^ 



4(^70 -m) (4^2 -m^^ 



(4^70 - m) + f^- ^ (1 - I sin(e/2)|) 

\Pi + Pfr 



MPi) 

(B2) 

which yields a second-order contribution to the cross section given by p7|) . As in the 
Coulomb case, the imaginary part of the scattering amplitude contains a divergent term, 
but the resulting polarization is finite at first-order. 

The second order scattering amplitude, A^2 = ^s(P/)f 2^tr(Pi)5 is defined by 

'^'^lBip,,k)-l±^Bik,p,), (B3) 



(27r)^ ■' kP' — w? + ie 

where 

B{P2,Pi) = ^ 4^70- ^1- ^2 - -^-^{P2 -Pi 

+ [iP^' - Pi'ho - 2E{^2- ^1)] , (B4) 

and q = P2 — Pi- Our task is made easier by splitting the calculation into several parts, so 
that 

B{P2,Pl) = Ciri(p2,Pl) +C2r2(P2>Pl) +^C3f3(p2,Pl), (B5) 



where 



,1) ci = -2ttGM, Ti{p2,Pi) = ^ 



(2) C2 = -47rGM, T2(p2,Pi) 



\Q\ 

iP2' - Pl'){^2- ^1] 

\q\' 



(3) tC3=i7rGM, T^ijP2,Pi) = ^ • (B6) 



We now consider each pair of terms in ()B3|) in turn, and wherever possible employ the 
simplification that ^iUriPi) = mUr{Pi) ^sijPf) = Us{pf)m. 

First, let us consider term (2x2), to illustrate a method for evaluating these integrals. 

^ ^ " ^ y (2vr)3 \k-p,\^\k-p,\^- ^^^^ 
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We start by defining Q = (py. — Pj)/2, R = (pj + py)/2 and moving to the centre of mass 
frame, where the 1 axis is along R and the 3 axis is along Q. Then 

fc I— > fe + R, 

fe^ - ^ fc^ - Q2 + 2R • k. (B8) 

Next, we employ spheroidal coordinates, {u,v,(l)}, so 

ki — |Q| sinhwsini>cos0, 
^2 = I Q I sinh u sin v sin 0, 

^3 = |Q| coshti cosv. (B9) 

With these coordinates the measure of integration is 

d^k = iQl^sinh u sin V (sinh^ u + sin^ v) dudydcf) (BIO) 

and the important quantities in the integral become 

\k — P2\\k — Pi\ ^ \k — Q\\k + Q\ = |Qp(sinh^ti + sin^t') 

k'^ -Q'^ ^\Q\'^{smh'^u-sm^v). (Bll) 

With these replacements, ^ — m ^ ft — Ejq — m —Jt — Ejq. This has only a spatial 
component, so will couple to the odd part R • fe of the rest of the integral. Therefore 

d'^k 4(fe'-Q2)(R.fe)(/(f-£;7o) 



(27r)3 |fe-Q|4|fc + Q|4 



C2^(m — £'7o) f , , sinh^usin^i'fsinh^u — sin^i)) /t^^^n 
' dudv — — 2"-^^ „ -. (B12) 



2|Q|7r2 J (sinh^ii + sin^f)^ 

The integral evaluates to 7r^/16, and the result is 

2f<2nf2 

(2x2) = ^— -^(m-£;7o). (B13) 
Similar techniques can be used to evaluate terms (3x3), and (1 x 2) + (2 x 1). We find 

(3x3) = ^'^'^\ {m - 7^70), (B14) 
(1 X 2) + (2 X 1) = — -^(9^70 - 5m). (B15) 
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The other cross terms, (1 x 3) and (2x3), can be shown to be zero by symmetry consid- 
erations. The remaining term, (1 x 1), has two parts: an integral which is tractable using 
the above techniques, and an integral with a pole at |fc| = |p| that is familiar from the 



second-order Coulomb calculation of Dalitz 

"(m-2E7o) 



;i X 1) = {27r'G'M^) 



\Pf-P. 



2l|. Explicitly, 

+ (2E7o-m)=^(m + E7o)(//7r3) 



(B16) 



where / and J are defined in terms of a low-frequency cut-off A as 



J 



2 y [(fc-p^)2 + A2][(fc-p.)2 + p](^2_ ^2+,g)- 

To lowest order in A these integrals evaluate to (see Itzykson and Zuber |41|]) 



(B17) 



/ 2p|sin(g/2)| \ 



J = sec^(^/2)/ + 



V cos2(0/2) ~ csc(0/2)) - cos2(0/2) 



ln(A/2p) 



I -J 



TT 



(csc(^/2) - 1) + 



■ln(sin2(0/2)). 



(B18) 



Vcos2(^/2)' ' ' Vcos2(^/2) 

The imaginary parts of / and J diverge as A ^ 0. However we note that the difference I — J 
is finite. We can rewrite the J and J parts of (1 x 1) as 



2/(2^2 - m^){2E-fo - m) + (/ - J){4:E^ - m^){m - E70). 



(B19) 



The 2i?7o — m associated with the / term ensures that this term does not contribute to the 
polarization at first order, as we see below. 

The result of summing (IbT3|) . (IbTI) . (IbT5|1 and (IbT6|i is 

n^G^M^l,.^ . 4(^70 -m) (4^2 -m2) 



r, = - 



(4^70 - m) + 



sin(0/2)|) 



ttG^M^AE^ - m^) ln(sin^(^/2)) 



m — Ejq) + ioo{2E'yQ — m). 



2p^ cos2(^/2) 

The second order contribution to the scattering cross section is 



(B20) 



dnj. 



27r/ 2 



2m 



2m 



m 



2m 



2m 



(B21) 
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where F = 7of ''"70. After employing the appropriate gamma matrix trace theorems, this 
evaluates to (fTTj) . The imaginary part of ()B20|) does not contribute to the total unpolarized 
scattering cross section, but it does contribute to the net polarization. The polarized cross 
section in direction s is given by 

^) = (^) ' 1^ {-"^ ^' ^ + ^ "-^^ ' 

= —I ^ - Im / - J)s ■ X p^). B22 

Here we have made use of the result 

Tr {75 ^ ^ (/ 4 = -Aieap^sa'^b'^c'd^, (B23) 

where ea/s-ys is antisymmetric under exchange of any pair of indices. The fraction of polarized 
flux V is found by dividing polarized flux ()B22|) by the unpolarized flux at first order ()16|) 
to obtain result (fTH|). 
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